Real-time gauge theory simulations from stochastic 
quantization with optimized updating 
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Stochastic quantisation is applied to the problem of calculating real-time evolution on a Min- 
kowskian space-time lattice. We employ optimized updating using reweighting, or gauge fixing, 
respectively. These procedures do not affect the underlying theory, but strongly improve the 
stability properties of the stochastic dynamics. 
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1. Introduction 

First-principles simulation of quantum field theories (such as QCD) is a notoriously 
hard problem of theoretical physics. Lattice calculations tipically use euclidean formula- 
tion, where one can apply importance sampling. In contrast, in Minkowski space, this is 
ineffective since the probability weight is in this case a phase factor: exp(iS'). Stochas- 
tic quantisation pQ, however can be generalised to complex weights. Recently stochastic 
quantisation was applied to calculate the real-time evolution of quantum fields f [2l3| and 
references therein). Recently stochastic quantisation has also been applied to nonzero 
chemical potential problems ^ . In [5] , optimized updatings were studied, which improve 
the behaviour of the complex Langevin equations, whitout changing the underlying the- 
ory. The main insight is the usage of the fixedpoints of the Langevin flow as a criteria 
for convergence. 

2. Stochastic Quantisation 

The time evolution can be formulated using the path integral formalism as an average 
weighted with exp(iS'): 
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Using stochastic quantization the real-time quantum configurations in 3+1 dimensions 
are constructed by a stochastic process in an additional (5th) Langevin-time 6, by use of 
a Langevin equation: 

d<^{9) SS 

where rj, the noise term satisfies: (?7(^)) — 0, {ri{9)rj{0' )) = 2S{9~9'). The expectation 
value of any observable can be calculated from the Langevin-time evolution using the 
following formula: 

T 

{om = ^J o{m)de (3) 



We applied stochastic quantisation to "one plaquette" models, which serve to illustrate 
and test concepts of optimized updating, a scalar oscillator (0-1-1 dimensional scalar 
field theory), pure 5/7(2) gauge theory, using the Wilson action on a 3+1 dimensional 
lattice. The field theories were discretised on a comlex time contour, which enabled us to 
either calculate equilibrium distributions, or non-equilibrium time evolution. The contour 
tipically has a downward slope, which improves convergence. 

The complexity of the drift term in ^ means that the originally real fields are com- 
plexified. This transfroms a real scalar field to complex scalar, the link variables of an 
SU{N) field theory to SL{N, C). Only after taking noise or Langevin-time averages, re- 
spectively, the expectation values of the original real scalar theory {SU (N) gauge theory) 
are to be recovered. Accordingly, if the Langevin flow converges to a fixed point solution 
of Eq. (21) it automatically fulfills the infinite hierarchy of Dyson-Schwinger identities of 
the original theory. 

3. Optimized updating of toy models 

As an example to illustrate reweighting we consider the one-plaquette model with U{1) 
symmetry. For U = e^"^ the action is given by 

5*0 = ^{U + U-^) = Pcos^ (4) 

with real coupling parameter (3. This simple model enables us to check the results of the 
complex Langevin equation by evaluating the integral 

(O)o = yJ ^^^^ ■ 



For real /3 the integrand in Eq. Q is not positive definite, which mimics certain aspects 
of more complicated theories in Minkowskian space-time. Solving the discretised Langevin 
equation corresponding to this model, using (3=1, and a Langevin stepsize e = 10~^, 
one gets: 



vithout 



e*'^)o -0.009(±0.006)-i0.00006(±0.00007). (6) 
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Evaluating the numerical integral in ([5]), one gets (e"^) ~ iO.575. One observes that 
the simulation yields a wrong result that is compatible with zero, in contrast to the 
non- vanishing imaginary value obtained analytically. 

The same averages may be calculated using by performing reweighting: changing the 
action, and recompensating the change in the measurable so that the end-result remains 
the same. 

Jg d(p e"-^" eH'So-s„) 

Evaluating the denominator and numerator of this factor one considers Sa as the action 
instead of 5o, thus the Langevin process is changed [5j. We will consider here the family 
of actions Sa = Sq + aip ~ (3 cos Lp + aLp. One finds that this new class of actions gives 
exact results for some region of its parameters (see Fig.[T]), such that (e*''')o is recoverable, 
using jll). 
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Fig. 1. The real and imaginary part of the average (e"''}a=i as a function of 0. The Hnes represent averages 
obtained from direct integration, while the symbols are measurements using a stochastic process. 

Qualitative understanding of the behavior of stochastic processes corresponding to 
Sa is possible by studying the fixed point structure of the drift term of the Langevin 
equation. In Fig. [2l the flowchart (normalized drift vectors) is plotted alongside a scatter 
plot of Lp on the complex plane. One observes (comparing with Fig.[T]) that the stochastic 
results brake down when the fixed point structure changes. Exact results correspond to 
an attractive fixedpoint, and a relatively compact distribution of field values, while the 
absence of attractive fixedpoints, and a wide distribution corresponds to the breakdown. 
The breakdown of the process using the original Sq action is consistent with this picture: 
it has no attractive fixedpoints. 

Consider now the SU{2) one plaquette model, as a second example. The action is given 
by: 

S{U) = ^TtU, (8) 

which is invariant under the "gauge" transformation U — > W~^UW with U,W G 
SL{2,C) (after complexification). 
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Fig. 2. Shown is the real and imaginary part of dSa/dip plotted as a vector with origin at each 99- value 
and with normalized length for better visibility, using a = 1, and /3 = 0.5 on the left panel and /3 = 1.5 
on the right panel. Also shown is the distribution of f as obtained from the full solution of the respective 
complex Langevin equation. 
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Fig. 3. The real and imaginar'y'paJt of the plaquette variable TrC/ from snapsliofs Hlfith constant Langevin— 
time stepping for the SU{2) one-plaquette model with 13 = 1. The left graph shows the wide distribution 
of values obtained from the standard Langevin dynamics, while the right graph displays the compact 
distribution from the (" gauge-fixed" ) optimized process (see text for explanation). 

Since the SU{2) one-plaquette model has a global symmetry, one may use this symme- 
try to "gauge-fix" certain variables in order to constrain the growth of fluctuations. In 
the following we will use it in order to diagonalize U after each successive Langevin-time 
step: U' = diag (a + iVi — a'^,a — iVl — a^) (for details, see [S]). The right graph of 
Fig. [3] shows that the "gauge-fixing" leads to a compact distribution, and exact results, 
in contrast to the distribution from the unmodified process displayed on the left of that 
figure, which gives wrong results. 

In conclusion, we have demonstrated the usage of optimized updating on specific ex- 
amples, which change the Langevin process by reweighting or using the symmetries of 
the theory, in a way that it gives exact results for the original theory. 
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